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Abstract. Factorial Schur functions are generalizations of Schur functions that have, in addition to 
the usual variables, a second family of "shift" parameters. We show that a factorial Schur function 

i— i times a deformation of the Weyl denominator may be expressed as the partition function of a 

particular statistical-mechanical system (six vertex model) . The proof is based on the Yang-Baxter 

u equation. There is a deformation parameter t which may be specialized in different ways. If t — — 1, 

then we recover the expression of the factorial Schur function as a ratio of alternating polynomials. 
If t = 0, we recover the description as a sum over tableaux. If t — oo we recover a description of 
Lascoux that was previously considered by McNamara. We also are able to prove using the Yang- 
Baxter equation the asymptotic symmetry of the factorial Schur functions in the shift parameters. 
Finally, we give a proof using our methods of the dual Cauchy identity for factorial Schur functions. 

|- Thus using our methods we are able to give thematic proofs of many of the properties of factorial 

OO Schur functions. 

o 

CO 

1 Introduction 

Factorial Schur functions are generalizations of ordinary Schur functions s\(z) = 
s\(zi, - ■ ■ , z n ) for which a surprising amount of the classical theory remains valid. In 
addition to the spectral parameters z — (zi, ■ ■ ■ , z n ) and the partition A they involve 
a set a = (oti, a 2 , a 3 , • • • ) of shifts that can be arbitrary complex numbers, and are 
denoted s\(z\a). In the original discovery of Biedenharn and Louck [2j, only the 
special case where a n = 1 — n was considered. Their original motivation, inspired by 
questions from mathematical physics, was to the decomposition of tensor products 
of representations with using particular bases. 

Other early foundational papers are Chen and Louck [5], who gave new founda- 
tions based on divided difference operators, and Goulden and Hamel [7] where the 



1 



analogy between Schur functions and factorial Schur functions was further developed. 
In particular they gave a Jacobi-Trudi identity See Louck [15] for further historical 
remarks. 

Biedenharn and Louck (in the special case a n = 1 — n) defined s\(z\a) to be a 
sum over Gelfand-Tsetlin patterns, and this definition extends to the genral case. 
Translated into the equivalent language of tableaux, their definition is equivalent to 
([7]) below. It was noticed independently by Macdonald [32] and by Goulden and 
Greene [6] that one could generalize the factorial Schur functions of Biedenharn 
and Louck by making use of an arbitrary set a of shifts. Macdonald observed an 
alternative definition of the factorial Schur functions as a ratio of two alternating 
polynomials, generalizing the Weyl character formula. This definition is ^ below. 

Both Macdonald and Goulden and Green also noticed a relationship with what 
are called supersymmetric Schur functions. These are Schur functions in two sets 
of variables, z = (21,22, •••) and w = (wi,Wz, ••■). They are called Hook Schur 
functions by Berele and Regev pQ and Schur- Littlewood functions by Bump and 
Gamburd [3]. They are defined in terms of the ordinary Schur functions by 



where c^ v is the Littlewood- Richardson coefficient, /x and v run through partitions 
and v' is the conjugate partition. The relationship between the factorial Schur func- 
tions and the supersymmetric Schur functions is this: although the s\(z\a) are sym- 
metric in the Zi they are not symmetric in the a.{. Nevertheless, as the number n of 
the spectral parameters Z{ tends to infinity, they become symmetric in the oti in a 
certain precise sense, and in the limit, they stabilize. Thus in a suitable sense 



Another important variant of the factorial Schur functions are the shifted Schur 
functions that were proposed by Olshanskii, and developed by Okounkov and 01- 
shanskii [23], [22]. Denoted s* x (xi,--- ,x n ), they are essentially the same as the 
factorial Schur functions of Biedenharn and Louck, but incorporate shifts in the pa- 
rameters so that they are no longer symmetric in the usual sense, but at least satisfy 
the stability property s x (xi, • • • , x n ) = s x (xi, • • • , x n , 0). These were applied to the 
representation theory of the infinite symmetric group. 

Molev and Sagan [20] give various useful results for factorial Schur functions, 
including a Littlewood-Richardson rule. A further Littlewood-Richardson rule was 
found by Kreiman [UJ. Knutson and Tao [TU] show that factorial Schur functions 




lim s\(z\a) = s\(z\\a). 
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correspond to Schubert classes in the equivariant cohomology of Grassmanians. See 
also Mihalcea [19] and Ikeda and Naruse [9]. 

Tokuyama [21] gave a formula for Schur functions that depends on a parameter 
t. This formula may be regarded as a deformation of the Weyl character formula. 
It was shown by Hamel and King [8] that Tokuyama' s formula could be generalized 
and reformulated as the evaluation of the partition function for a statistical system 
based on the six-vertex model in the free-fermionic regime. Brubaker, Bump and 
Friedberg [3] gave further generalizations of the results of Hamel and King, with new 
proofs based on the Yang-Baxter equation. More specifically, they showed that the 
six-vertex model in the free-fermionic regime satisfies a parametrized Yang-Baxter 
equation with nonabelian parameter group GL(2) x GL(1). Using this, they showed 
how to give statistical-mechanical systems whose partition functions were Schur func- 
tions times a deformation of the Weyl denominator. This result generalizes the results 
of Tokuyama and of Hamel and King. 

Lascoux [12], [13] found six- vertex model representations of Grassmannian Schu- 
bert polynomials. A proof of this representation based on the Yang-Baxter equation 
was subsequently found by McNamara [IS] . 

We will give a simultaneous generalization of [3] and of the representations of Las- 
coux and of McNamara. As in [3] we will consider partition functions of statistical- 
mechanical systems in the free-fermionic regime. A significant difference between 
this paper and that was that in [3] the Boltzmann weights were constant along the 
rows, depending mainly on the choice of a spectral parameter z^. Now we will con- 
sider systems in which we assign a spectral parameter Zi to each row, but also a 
shift parameter ct,- to each column. Furthermore, we will make use of a deformation 
parameter t that applies to the entire system. 

Using the Yang-Baxter equation, we will show that the partition function may be 
expressed as the product of a factor (depending on t) that may be recognized as a de- 
formation of the Weyl denominator, times the factorial Schur function s\(z\a). This 
is thus a Tokuyama-type formula for the factorial Schur functions. By specializing 
t in various ways, we will obtain different formulas for the factorial Schur functions. 
Taking t = 0, we obtain the tableau definition. Taking t — —1, we obtain the repre- 
sentation as a ratio of alternating polynomials, which was Macdonald's generalization 
of the Weyl character formula. Next, when t = oo we recover the representation of 
Lascoux and McNamara. Although we do not prove the supersymmetric limit (JIJ, 
we will at least prove the key fact that the s\(z\at) are asymptotically symmetric in 
the ctj as the number n of spectral parameters Zi tends to infinity. We will obtain 
this by another application of the Yang-Baxter equation. We also give a proof of the 
dual Cauchy identity for factorial Schur functions using our methods. 
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2 Factorial Schur functions 

A good primary reference for factorial Schur functions is Macdonald [TB]. They are 
also in Macdonald [T7j, Ex. 20 in Section 1.3 on p. 54. 

Let oti, Q!2, Q!3, • • • be a sequence of complex numbers. If z £ C let 

= (z + ai) • • • (z + a r ). 

Macdonald [16J gives two formulas for factorial Schur functions that we will also 
prove to be equivalent by our methods. Let /i = (fix,-- - , fi n ) where the /i, are 
nonnegative integers. Let z±, ■ ■ ■ , z n be nonzero complex numbers which we will call 
spectral parameters. Define 

A^(z\a) = det((zi\a)^) i>:j 
where 1 ^ i, j ^ n in the determinant. We will also use the notation 



Let 8 — (n — 1, n — 2, ■ ■ ■ ,0) and let A = (Ai, • • ■ , A n ) be a partition of length at 
most n. Define 

/ i \ A x+ s{z\a) 
s\{z\a) = . 2 

A s {z\a) 

The denominator here is actually independent of a and is given by the Weyl denom- 
inator formula: 

A&{z\a) = Y[(zi - Zj) (3) 

i<j 

Indeed, it is an alternating polynomial in the Zi of the same degree as the right-hand 
side, and so the ratio is a polynomial in a that is independent of z*. To see that the 
ratio is independent of a, one may compare the coefficients of z s on both sides of 

©■ 

Theorem 1 The ratio s\(z\a) is a symmetric polynomial in zx,- ■ ■ , z n . 
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We will prove this at the end of Section [3] 

Macdonald also gives a combinatorial formula for the factorial Schur function as 
a sum over semi-standard Young tableaux of shape A in {1, 2, ■ • • , n}. Let T be such 
a tableau. If i, j are given such that j ^ Aj, let T(i,j) be the entry in the i-th row 
and j-th column of T. Let 



T*(i,j) = T(i,j)+j-i, 



and define 



(z\a) T = Y[(z T (i,j) + aT*(ij)). 

(m) 

For example, if A = (4, 2, 0) and n — 3, then we might have 



(4) 
(5) 
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(6) 



and 



(z\a) T = (zi + a l )(zi + a 2 )(zi + a 3 )(z 3 + a 6 )(z 2 + a 1 )(z 2 + a 2 ). 



Theorem 2 Let X be a partition. Then 

s x (z\a) = ^2(z\a) T } (7) 

T 

where the sum is over semistandard Young tableaux with shape A in 1, 2, 3, • • • , n. 

We will prove this in Section |4j If A is a partition then A' will denote the conjugate 
partition whose Young diagram is the transpose of that of A. 
Let /i be a partition, and let a M denote the sequence 



Theorem 3 (Vanishing Theorem) We have 



S\[-Oiu \Oi) 



if A £ //, 



See Okounkov [2T] Section 2.4 and Molev and Sagan |20j. We will prove this in 
Section [U 
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3 Yang-Baxter equation 



We review the six- vertex model and a case of the Yang-Baxter equation from [3j. 
We will consider a planar graph. Each vertex v is assumed to have exactly four 
edges adjacent to it. Interior edges adjoin two vertices, and around the boundary 
of the graph we allow boundary edges that adjoin only a single vertex. Every vertex 
has six numbers ai(v), 02(1;), 61 (v), b2(v), C\(v), 02(f) assigned to it. These are called 
the Boltzmann weights at the vertex. Once we have specified the graph, the Boltz- 
man weights at the vertices, and the boundary spins, we have specified a statistical 
system &. 

A state s of the system will be an assignment of spins to the interior edges. Given 
a state of the system, every edge, boundary or interior, has a spin assigned to it. 
Then every vertex will have a definite configuration of spins on its four adjacent 
edges, and we assume these to be in one of the two orientations listed in (J8J). Then 
let s (v) equal ai(v ), a 2 {y ), &i(t> ), b 2 {y ), C\{v ) or c 2 {v ) depending on the configuration 
of spins on the adjacent edges. If v does not appear in the table, the weight is zero. 
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The Boltzmann weight of the state /3(s) is the product Y[ v Ps{ v ) of the Boltzmann 
weights at every vertex. We only need to consider configurations in which the spins 
adjacent to each vertex are in one of the configurations from the above table; if this 
is true, the state is called admissible. A state that is not admissible has Boltzmann 
weight zero. 

The partition function Z(&) is ^2 s f3{s), the sum of the Boltzmann weights of 
the states. We may either include or exclude the inadmissible states from this sum, 
since they have Boltzmann weight zero. 

If at the vertex v we have 

ax{v)a 2 (y) + b x (v)b 2 {v) - Ci(v)c 2 (v) = 0, 

the vertex is called free-fermionic. We will only consider systems that are free- 
fermionic at every vertex. 
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Brubaker, Bump and Friedberg proved the following version of the Yang-Baxter 
equation. 

Theorem 4 Theorem 3] Let v, w be vertices with free-fermionic Boltzmann 
weights. Define another type of vertex u with 

ai(u) = a 1 (v)a 2 (w) + b 2 (v)b 1 (w), 

a 2 (u) = bi(v)b 2 {w) + a 2 (v)a 1 (w) , 

bi(u) = b 1 {v)a 2 (w) - a 2 (v )bi(w), 

b 2 (u) = —ai(v)b 2 (w) + b 2 (v)a 1 (w), 

Ci(«) = ci(w)c 2 (w), 

c 2 {u) = c 2 {v)ci{w). 

Then for any assignment of edge spins £j G {±} (i = 1, 2, 3, 4, 5, 6) the following two 
configurations have the same partition function: 




Note that by the definition of the partition function, the interior edge spins 
(labeled v, 7, \i and 5, ip, 4>) are summed over, while the boundary edge spins, labeled 
Si are invariant. In order to obtain this from Theorem 3 of j3] one replaces the 
R-matrix 7r(R) in the notation of that paper by a constant multiple. 

Now let us consider two types of Boltzmann weights. Let Z\ 1 ■ ■ ■ ,z n be spec- 
tral parameters, oti, a 2 , a^, ■ ■ ■ another sequence of complex numbers, and t another 
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parameter. If 1 ^ i, k ^ n and if j ^ are integers, we will use the following weights. 
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tZi + Zfc 



Z k (t + l)Zi 



(t + l)z k 



Lemma 1 We may take u = vrr{h k,t) , v = Vr(i,j,t) and w = Vr(k,j,t) in Theo- 
rem^ 

Proof The relation a^u) = ai(v)a 2 (w) + b 2 (v)bi(w) becomes 

tZi + z k = 1 ■ (z k - taij) + {zi + aj)t, 

and all the other relations are checked the same way. □ 
We will make use of two special n-tuples of integers, namely 

p= (tv ,3,2,1), <J=(n-l,n-2,... ,2,1,0). 

Let A = (Ai, • • • , A n ) be a partition, so Ai ^ . . . ^ A n ^ 0. Let us consider a lattice 
with n rows and n + X\ columns. We will index the rows from 1 to n and the columns 
from 1 to n+ Ai, in reverse order. We put the vertex vr(i,j,t) at the vertex in the 
i row and j column. 

We impose the following boundary edge spins. On the left and bottom boundaries, 
every edge is labeled +. On the right boundary, every edge is labeled — . On the 
top, we label the edges indexed by elements of form X 3 ; + n — j + 1 for j = 1, 2, • • • , n 
with a — sign. These are the entries in A + p. The remaining columns we label with 
a + sign. 

For example, suppose that n = 3 and that A = (5, 4, 1), so A + p = (8, 6, 2). Since 
Xj + n — j + 1 has the values 8, 6 and 2, we put — in these columns. We label the 
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vertex vr(i,j,t) in the i row and j column by ij as in the following diagram: 
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<±> (10) 



Let this system be called ©£ t . 

We recall that a Gelfand-Tsetlin pattern is an array 



1 = < 



Pll Pl2 
P22 



Plr. 



P2n 



(11) 



Pr, 



in which the rows are interleaving partitions. The pattern is strict if each row is 
strongly dominant, meaning that pa > Pi^+i > ■ ■ ■ > Pi n . 

A staircase is a semistandard Young tableau of shape (Ai + n, X± + n — 1, . . . , Ai)' 
filled with numbers from {1, 2, ... , Ai +n} with the additional condition that the di- 
agonals are weakly decreasing in the south-east direction when written in the French 
notation. 

An example of a staircase for A = (5,4, 1) is 
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Proposition 1 Let X be a partition. There are natural bisections between the fol- 
lowing three sets of combinatorial objects 

1. States of the six vertex model &\ t - 

2. Strict Gelfand-Tsetlin patterns with top row A + p, and 

3. Staircases whose rightmost column consists of the integers between 1 and Ai + n 
that are not in A + p. 

Proof Suppose we start with a state of the six vertex model <5^ t . Let us record 
the locations of all minus signs that live on vertical edges. Between rows k and 
k + 1, there are exactly n — k such minus signs for each k. Placing the column 
numbers of the locations of these minus signs into a triangular array gives a strict 
Gelfand-Tsetlin pattern with top row A + p. 

Given a strict Gelfand-Tsetlin pattern 1 = (t^) with top row A + p, we construct 
a staircase whose rightmost column is missing A + p in the following manner: We fill 
column j + 1 with integers ui, . . . , u n ^j + \ 1 such that 



We give an example. As before let n — 3 and A = (5, 4, 1). Here is an admissible 
state: 




Then the entries in the i-th row of X are to be the columns j in which a — appears 
above the vertex. Therefore 



{1, 2, . . . , n + Ai} — {ui, . . . , Wn-j+Ai} U {t n +i-j,i, ■ ■ ■ t n +i-j,j-i}- 



It is easily checked that these maps give the desired bijections. 
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Taking the complements of the rows (including a fourth empty row) gives the fol- 
lowing sets of numbers: 

7 5 4 3 1 

8 6 5 3 2 1 

8 7 6 5 3 2 1 
8 7 6 5 4 3 2 1 



so the corresponding staircase is (12) 
Proposition 2 The function 



Y[(tZj + Zi) 



.i>3 



(13) 



is symmetric under permutations of the Zi 



Proof Let & 



It is sufficient to show that (13) is invariant under the 



interchange of z^ and z i+ i. The factors in front are permuted by this interchange 
with one exception, which is that tzi + z i+ i is turned into tz i+ i + Zj. Therefore we 
want to show that (tzi + Zi + \)Z{&) = {tz i+ i + Zj)Z{&) where & is the system 
obtained from & by interchanging the weights in the i, i + 1 rows. Now consider the 
modified system obtained from & by attaching vrr(i, i + l,t) to the left of the i, i + 1 
rows. For example, if i — 1, this results in the following system: 
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Referring to ([8]) , there is only one admissible configuration for the two interior edges 
adjoining the new vertex, namely both must be +, and so the Boltzmann weight of 
this vertex will be tz^ + z i+ i. Thus the partition function of this new system equals 
(tzi + z i+ i)Z(&). Applying the Yang-Baxter equation repeatedly, this equals the 
partition function of the system obtained from & by adding vrr(i,i + to the 
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right of the i, i + 1 rows, that is, {tz i+1 + Zi)Z(&), as required. This proves that (13) 
is symmetric. □ 



Theorem 5 We have 



Proof We show that the ratio 



\\itZj + Zi 



s x (z\a). 



UiKjitZj +Zi) 



(14) 



is a polynomial in the Z{, and that it is independent of t. Observe that (13) is 
an element of the polynomial ring C[zi, ■ ■ ■ ,z n ,t], which is a unique factorization 
domain. It is clearly divisible by tzj + Zi when i > j, and since it is symmetric, it 
is therefore divisible by all tzj + zi with i ^ j. These are coprime, and therefore it 



is divisible by their product, in other words (14) is a polynomial. We note that the 



numerator and the denominator have the same degree in t, namely ^n(n — 1). For 
the denominator this is clear and for the numerator, we note that each term is a 
monomial whose degree is the number of vertices with a — sign on the vertical edge 
below. This is the number of vertical edges labeled — excluding those at the top, that 
is, the number of entries in the Gelfand-Tsetlin pattern in Lemma [T] excluding the 
first row. Therefore each term in the sum Z{& T X t ) is a monomial of degree \n{n — 1) 
and so the sum has at most this degree; therefore the ratio (14) is a polynomial of 
degree in t, that is, independent of t. 

To evaluate it, we may choose t at will. We take t = —1. 

Let us show that if the pattern occurs with nonzero Boltzmann weight, then 
every row of the pattern (except the top row) is obtained from the row above it by 
discarding one element. Let /i and v be two partitions that occur as consecutive rows 
in this Gelfand-Tsetlin pattern: 



/'2 



V2 



where /i^ are the column numbers of the vertices (i,Hk) that have a — sign on the 
edge above the vertex, and are the column numbers of the vertices (i, v^) that 
have a — sign on the edge below it. Because t — — 1, the pattern 



e- 
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does not occur, or else the Boltzmann weight is zero, and the term may be discarded. 
Therefore every — sign below the vertex in the i-th row must be matched with a — 
above the vertex. It follows that every z/j equals either /Xj or //j+i. Thus the partition 
v is obtained from \i by discarding one element. 

Let Hk be the element of \i that is not in v. It is easy to see that in the horizontal 
edges in the i-th row, we have a — sign to the right of the /x^-th column and + to 
the left. Since t = —1, the Boltzmann weights of the two patterns: 



e- 



e 







e- 





both have the same Boltzmann weight Zi + atj. We have one such contribution from 
every column to the right of /ik and these contribute 

Hizi + aj) = tela)"*" 1 . 
j'=i 

For each column j = /ii with I < k we have a pattern 



and these contribute —1. Therefore the product of the Boltzmann weights for this 
row is 

(-l^tela)"*" 1 . 

Since between the i-th row /i and the (i + l)-st row v of the Gelfand-Testlin 
pattern one element is discarded, there is some permutation a of {1,2,3,- •• , n) 
such that v is obtained by dropping the <r(i)-th element of /i. In other words, 

fi k ~ 1 = (A + - 1 = (A + 5)^(1) 

and we conclude that 

^A,-i)=E ± n^i«) (A+5)CTW - 
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The sign may be determined by taking all aj = 0. The term corresponding to a is 



then ±n i ^ + ^ <rW . Because the ratio (14) is symmetric, and since (when t = — 1) 



the denominator is antisymmetric, we have 



Z{&l_ x ) = ±J2(- 1 ) li ' T) W {X+pUi) = ±Ax+s(z\a). 



<t&S„ 



Thus using ^ we see that (14) equals ±s\(z\a). To evaluate the sign, we may take 
t = and all z% — 1. Then the partition function is a sum of positive terms, and 
S\(l, ■ • • , 1) is positive, proving that (14) equals s\(z\a). □ 



Proof of Theorem [T] Theorem [T] follows from Proposition [2] and Theorem [5j □ 



4 Proof of Theorem 2 



When t = 0, the Boltzmann weig htfor f is zero, so we are limited to states 
omitting this configuration. If T = Ta+ p is the Gelfand-Tsetlin pattern corresponding 



to this state, and if the entries of % in are denoted pn- as in (11), it is easy to see 
that the equality Pi-1^-1 = Pi,k would cause this configuration to appear at the i, j 
position, where j = p i; k- Therefore Pi-i^-i > Pi,k- This inequality implies that 
we may obtain another Gelfand-Tsetlin pattern Ta with top row A by subtracting 



Pn—i+l 



in 



i + 1, n 



1) from the i-th row of %\ +p . 



Consider the example of A = (4, 2, 0) with n = 3. Then 



if 1 



A+p 





We associate with Ta a tableau T(%\) of shape A. In this tableau, removing all 



boxes labeled n from the diagram produces a tableau whose shape is the second 
row of Ta- Then removing boxes labeled n — 1 produces a tableau whose shape is 



the third row of Ta, and so forth. Thus in the example, T(T\) is the tableau T in 
Let wq denote the long element of the Weyl group S n , which is the permutation 
i i— > n + 1 — i. 
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Proposition 3 Let t = 0, and let s = s(%\ +p ) be the state corresponding to a special 
Gelfand-Tsetlin pattern %\ +p . With Z\ as above and T = T(%\) we have 




The locations labeled o produce powers of z^ and the locations labeled • produce 
shifts of the form Zi + ctj. The weight of this state is 

n^s(^) = zfz 2 (z 1 + a 6 )(z 2 + a 2 ){z 2 + ai)(z 3 + a 3 ){z 3 + a 2 ){z 3 + 

Applying wq interchanges z-i < — > z±-i. The factor z\z 2 becomes z\z 2 = z w °( p \ and 
the terms that remain agree with 



Proof We are using the following Boltzmann weights. 
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We have contributions of Z{ from vertices that have — on the vertical edge below, and 
there is one of these for each entry in the Gelfand-Tsetlin pattern. These contribute 
a factor of z p . Applying wq to z leads to z Wo ^ 

Considering the contribution from the vertex, there will be a factor of zi + 

e. 

a.,- when the vertex has the configuration © —<H ^© . In the above example, the 



locations are labeled by •. 

Let Ta+p be the Gelfand-Tsetlin pattern (11). Let %\ and T = T(%\) be as 
described above, and let T* be as in Q. Let the entries in %\ +p and Ta be denoted 
Pi j and qij, with the indexing as in 11 Thus q^j = Pij — n + j — 1. We will make 
the convention that qi )Tl +i = Pi,n+i — 0. The condition for • in the i,j position of the 
state is that for some k with i ^ k ^ n we have Pi+x^+i < j < Pi,k- Translating this 
in terms of the qij = pij — n + j — 1 the condition becomes gj+i^+i + n — k + 1 ^ 
j < q%,k + n — k. The effect of wq is to interchange Zi -H- z n -i+i, and therefore 

n n 1i,k+n-k 




n n n (^+1-*+%) 

k=i j=q i+ i ik+1 +n-k+l 



(16) 



i=i 



n 



1 



On the other hand, in the tableau T, the location of the entries equal to 
in the (k + 1 — z)-th row is between columns q^+i^+i + 1 through q it k, and if j is one 
of these columns then 



T(k + l-j,j)=n + l-i, 
Therefore in the notation ^ we have 



T*(k + l-j,j)=n + j-k. 



Qi,k 



{z\a) rj 



n n n (z n+1 -i+a n+j - k ) 

k=i J=9i+i,fe+i+l 



i=l 



This equals (16) and the proof is complete. □ 

We now give the proof of Theorem [2j Summing over states, the last Proposition 
implies that 

Z x (w o (z),a,0) = 
Since s\(z,a, 0) = s\ (wq(z), a, 0) we have 



r Z A (w (z),a,0) s 



(*,a,0) 



'Wo(p) 



E 1 



2; a 



and the statement follows. 
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5 The limit as t tends to infinity 



In this section, we examine the behaviour of our six-vertex model as we send the 
parameter t to infinity. The first result of this section may be construed as a red- 
erivation of a theorem of Lascoux [T31 Theorem 1]. However the approach we take 
is to interpret it as giving us a proof of the equivalence of factorial Schur functions 
and double Schubert polynomials for Grassmannian permutations. We also obtain 
the vanishing theorem for factorial Schur functions (Theorem [3]). 

We start with two simple lemmas. To state these, we describe the six admiss- 
able arrangements of edge spins around a vertex as types a 1; a 2 , &i, b 2 , c\ and c 2 
respectively, when reading from left to right in the diagram ([8]). 

Lemma 2 In each state, the total number of sites of type a 2 , b\ and c\ is equal to 
n{n - l)/2. 

Proof This number is equal to the number of minus signs located in the interior of 
a vertical string. □ 

Let \x be the partiton (A + p)' . 

Lemma 3 In each state, the number of occurrences of a 2 , b 2 and C\ patterns in the 
i-th column is equal to /ij. 

Proof This number is equal to the number of minus signs located on a horizontal 
string between the strings labelled % and i + 1. This count is known since for any 
rectangle, knowing the boundary conditions on the top, bottom and rightmost sides 
determines the number of such spins along the leftmost edge. □ 

As a consequence of Lemma |2j the result of taking the limit as t — > oo can be 
interpreted as taking the leading degree term in t for each of the Boltzmann weights 
vr(i,j,t) as in Section [3] Thus with the set of Boltzmann weights (1, — ctj, 1, Zi + 
oij,Zi, 1) and corresponding partition function Z(&^ OQ ), we have 

Z(&U = z<>s x (z\a). 

By Lemma [3j if we consider our ice model with the series of Boltzmann weights 
in the following diagraml, then the corresponding partition function is given by 
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dividing Z{&\ ) by (-a) M . 



Gamma 
Ice 




B 






6 




6 


& 




5 


e- 












e- 




e- 


pe 


Boltzmann 
weight 


i 


i 


i 


-^i/ttj - 1 




i 



Let us denote the partition function for this set of Boltzmann weights by Z{<E)^ i00 {z\a)). 
Thus we obtain the result of [181 Theorem 1.1], 



Z(G 



A,oo 



[z\a)) 



7 s A (z|a). 



(17) 



We now pause to introduce the notions of double Schubert polynomials and Grass- 
mannian permutations so that we can make the connection to [13, Theorem 1] precise. 

A permutation is Grassmannian if it has a unique (right) descent. 

Let n and m be given and let A = (Ai, • • • , A n ) be a partition such that Ai ^ m. 
Then there is associated with A a Grassmannian permutation w\ G S n+m . This is 
the permutation such that 



K+i-i + i if i < n, 
i — X':„ if % > n, 



where A' is the conjugate partition. This has W\(n+1) < w\{n) and no other descent. 

Let Xi, ■ ■ ■ , x n+m and y%, ■ ■ ■ , y n+m be parameters. We define the divided differ- 
ence operators as follows. lil^i<n + m and / is a function of the Xi let 



d i f(x 1 , 



f-sj 



■I i 



where Sj/ is the function obtained by interchanging Xj and Xi + \. Then if w e S, : 



let w = Si t • • • Si k be a reduced expression of w as a product of simple reflections. 
Then let d w = ■ ■ ■ di k . This is well defined since the divided difference operators 
di satisfy the braid relations. Let wq be the long element of S n+m . Then the double 
Schubert polynomials, which were defined by Lascoux and Schutzenberger [H] are 
given by 



n & ~ 
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Theorem 6 The factorial Schur functions are equal to double Schubert polynomials 
for Grassmannian permutations. More precisely, 



&w x (x,y) = s x (x\ - y) 

Although this comparison result is well-known, it does not appear to have a 
good reference in the literature. One place where the relationship between factorial 
Schur functions and double Schubert polynomials for Grassmannian permutations is 
discussed is in Mihalcea 1191, Remark 2 in Section 5. 



Proof The proof is a comparison of (17) with [13, Theorem 1]. To translate the 



left-hand side in (17) into the staircase language of Lascoux we use the bijection in 
Proposition [TJ In Theorem 1 of [13] if Lascoux's x is our z and his y is our —a, then 
the left-hand side of his identity is exactly the partition function on the left-hand 



side of our (17). The monomial x pr y ^ on the right-hand side of his identity equals 



the monomial the right-hand side of fll7h. Lascoux' is the 

double Schubert polynomial & Wx . The statement follows. □ 

To conclude this section, we shall use this description to give a proof of the char- 



acteristic vanishing property of Schur functions. In lieu of (17) above, the following 
result is clearly equivalent to Theorem [3j 

Theorem 7 For two partitions A and \i of at most n parts, we have 

Z{G r X oQ {—a^\a)) = unless X C \i, 

^(®a!oo(-"a|«)) = 




Proof Fix a state, and assume that this state gives a non-zero contribution to the 
partition function Z{&\ 0O (— a^a)). Under the bijection between states of square ice 
and strict Gelfand-Tsetlin patterns, let ki be the rightmost entry in the z-th row of 
the corresponding Gelfand-Tsetlin pattern. We shall prove by descending induction 
on i the inequality 

n + 1 — % + Hi > ki. 

For any j such that ki+i < j < ki, there is a factor (xi/aij — 1) in the Boltzmann 
weight of this state. We have the inequality n+1 — i+f/* > n+1 — + fc+i > k i+ i 
by our inductive hypothesis. Since (a M )i = in order for this state to give 

a non-zero contribution to Z(&\ 00 (—a fl \a)), we must have that n + 1 — i + fa > k, 
as required. 



i ■ 
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Note that for all i, we have ki > n + 1 — i + \. Hence \ii > A; for all i, showing 
that fi D A as required, proving the first part of the theorem. 

To compute Z(©^ 00 (— a\\a)), notice that the above argument shows that there 
is only one state which gives a non-zero contribution to the sum. Under the bijection 
with Gelfand-Tsetlin patterns, this is the state with p t j = pij for all The 
formula for Z(©^ 00 (— a^lct)) is now immediate. □ 

6 Asymptotic Symmetry 

Macdonald [16j shows that the factorial Schur functions are asymptotically symmetric 
in the «j as the number of spectral parameters Zi increases. To formulate this prop- 
erty, let a be a permutation of the parameters a = a 2 , • • • ) such that <r(<x,) = <x,- 
for all but finitely many j. Then we will show that if the number of spectral param- 
eters Zi is sufficiently large, then s\(z\aa) = s\(z\a). How large n must be depend 
on both A and on the permutation a. 

We will give a proof of this symmetry property of factorial Schur functions us- 
ing the Yang-Baxter equation. In the following theorem we will use the following 
Boltzmann weights: 
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Theorem 8 Wii/i £/ie above Boltzmann weights, and with e±, ■ ■ ■ ,6$ fixed signs ±, 
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the following two systems have the same partition function. 





© © 



Proof This may be deduced from Theorem 3 of [3] with a little work. By rotating 
the diagram 90° and changing the signs of the horizontal edges, this becomes a case 
of that result. The R-matrix is not the matrix for u given above, but a constant 
multiple. We leave the details to the reader. □ 

Let Oi be the map on sequences a = (a±, «2, • • • ) that interchanges oii and 
Let A be a partition. We will show that sometimes: 



s x (z\a) = sx(z\aia) 



(18) 



and sometimes 



sx(z\a) = s x (z\o-ia) + s fi (z\a i a)(a i+1 - a { ), 



(19) 



where p is another partition. The next Proposition gives a precise statement distin- 



guishing between the two cases (18) and (19). 



Proposition 4 (i) Suppose that i + 1 G A + p but that i ^ A + p. Let p be the 

partition characterized by the condition that p + p is obtained from X + p by replacing 



the unique entry equal to i + 1 by i. Then (19) holds with this p. 



(ii) If either i + l^\ + poriE\ + p then (18) holds. 



For example, suppose that A = (3, 1) and n = 5. Then A + p = (8, 5, 3, 2, 1). If i = 4, 
then i + l = 5GA + p but i £ A + p, so (19) holds with p + p = (8, 4, 3, 2, 1), and so 
p=(3). 

Proof We will use Theorem [8] with a = a% and = oti + \. The parameter t may 
be arbitrary for the following argument. We first take i = n and attach the vertex 
u below the % and i + 1 columns, arriving at a configuration like this one in the case 
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n = 3, A 



(4,2,1). 
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There is only one legal configuration for the signs of the edges between u and the 
two edges above it, which in the example connect with (3,5) and (3,4): these must 
both be +. The Boltzmann weight at u in this configuration is unchanged, and the 
partition function of this system equals that of &\ t - After applying the Yang-Baxter 
equation, we arrive at a configuration with the u vertex above the top row, as follows: 
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26 
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24 



23 



22 



*H3 
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<( 38 (( 37 (| 36 (( 35 )( 34 )| 33 )( 32 ^^31 

©©(±)©©©©© 

Now if we are in case (i), the signs of the two edges are — , + then there are two legal 
configurations for the vertex, and separating the contribution these we obtain 



z(e x , t ) = (<TiZ(e x ,t)) + (0 - a) ((7^(6^)) 
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If we are in case (ii), there is only one legal configuration, so Z{&\,t) = ( a i^(^x,t)) ■ 

□ 

Corollary 1 Let X be a partition, and let I be the length of X. If n ^ I + i then 
s\(z\a) = s\(z\(Jia). 

Proof For if I is the length of the partition A, then the top edge in column j is a — 
sign when j ^ n — I, and therefore we are in case (i). □ 

This Corollary implies Macdonald's observation that the factorial Schur functions 
are asymptotically symmetric in the Oj. 

7 The Dual Cauchy Identity 

Let m and n be positive integers. For a partition A = (Ai, . . . , A n ) with Ai < m, we 
define a new partition A = (Ai, . . . , X m ) by 

Ai = \{j | Xj <m- i}\. 

After a reflection, it is possible to fit the Young diagrams of A and A into a rectangle. 
We shall prove the following identity, known as the dual Cauchy identity. 

Theorem 9 (Dual Cauchy Identity) For two finite alphabets of variables x = 
(a?!, . . . , x n ) andy = (y x ,. . . ,y m ), we have 

n m 

n + = yi s ^ x \ a ) s \(y\ - «)• 

The sum is over all partitions X with at most n parts and with X± < m. 

Proof The proof will consist of computing the parition funciton of a particular six 
vertex model in two different ways. We will use the weights vr(i,j,t) introduced in 
Section [3] with t = 1 and for our spectral parameters (zi, . . . , z m+n ), we will take the 
sequence (y m , ...,y x ,Xi,.. .,x n ). 

As for the size of the six-vertex model we shall use and the boundary conditions, 
we take a (m + n) x (m + n) square array, with positive spins on the left and lower 
edges, and negative spins on the upper and right edges. By Theorem [5] (with A = 0), 
the partition function of this array is 

+ + x ^ Tit** + 

i<j i,j i<j 
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We shall parition the set of all states according to the set of spins that occur 
between the rows with spectral parameters y\ and X\. Such an arrangement of spins 
corresponds to a partition A in the usual way, i.e. the negative spins are in the 
columns labelled Aj + n — i + 1. In this manner we can write our partition function 
as a sum 



E 



^top ^-bottom 
A A 



Here Z^ p is the parition funciton of the system with m rows, m + n columns, 
spectral parameters y m , . . . , y\ and boundary conditions of positive spins on the left, 
negative spins on the top and right, and the A boundary condition spins on the 
bottom. And ^bottom j g ^ e partition function of the system with n rows, m + n 
columns, spectral parameters xi, . . . ,x n and the usual boundary conditions for the 
partition A. 

By Theorem [5j we have 



7bottom 



\\{xi + Xj)s\(x\a). 



i<j 

It remains to identify Z^ p . To do this, we perform the following operation to 
the top part of our system. We flip all spins that lie on a vertical strand, and then 
reflect the system about a horizontal axis. As a consequence, we have changed the 
six-vertex system that produces the partition function Z^ p into a system with more 
familiar boundary conditions, namely with positive spins along the left and bottom, 
negative spins along the right hand side, and along the top row we have negative 
spins in those columns labelled Aj + m — % + 1. 

The Boltzmann weights for this transformed system are (recalling that t = 1) 
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Now again we use Theorem [5] to conclude that 

+ yj) s \(y\ 



— a). 



i<j 



Comparing the two expressions for the partition function of this square system com- 
pletes the proof. □ 
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